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ABSTRACT 

The recent detection bv lJee et all (|2009t ) of the massive cluster XMMU J2235. 3-2557 
at a redshift z « 1.4, with an estimated mass M324 = (6.4 ± 1.2) x 10^"* Mq, has 
been claimed to be a possible challenge to the standard ACDM cosmological model. 
More specifically, the probability to detect such a cluster has been estimated to be 
~ 0.005 if a ACDM model with gaussian initial conditions is assumed, resulting in 
a 3(7 discrepancy from the standard cosmological model. In this p aper we propose 
to use high redshift clusters as the one detected in iJee et"aLl (|2009D to compare the 
cosmological constant scenario with interacting dark energy models. We show that 
coupled dark energy models, where an interaction is present between dark energy and 
cold dark matter, can significantly enhance the probability to observe very massive 
clusters at high redshift. 
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1 INTRODUCTION 

One of the main challenges of present cosmology is to use 
observations to distinguish between a cosmological constant 
scenario and dynamical dark energy. More generally, it is 
a crucial task to devise new observational tests capable of 
detecting possible failures of the standard ACDM model, 
also in view of future observations, and to constrain the 
parameter space of alternative scenarios. The existence of 
massive clusters at high redshift has been recently used to 
test ACDM cosmologies. In pa rticular, weak lensing obser- 
vatio ns in a survey of 11 deg'^ (|jee et al.ll2009l : lRosati et al.l 
I2OO9I ) have recently detected a cluster of mass M324, = 
(6.4± 1.2) X IO^Mq (M200 = (7.3± 1.3) x IO^Mq) at a red- 
shift z ~ 1.4. The observation of such a massive cluster at 
high redshift represents an extremely rare event in the con- 
text of hierarchical structure formation. More specifically, 
the probability to detect such a cluster within a ACDM 
model with gaussian initial conditio ns in the v olume con- 
sidered by the sur vey presented in |jee et al] l|2009l 'l and 
iRosati et all l|2009l ) has been estimated by the same au- 
thors to be ~ 0.005, resulting in a 3a discrepancy with 
the assumed model. Naturally, this discrepancy can indi- 
cate that either initial conditions are not perfectly gaus- 
sian, or that the ACDM model itself has to be modified. 
Some analysis have already been carried out in order to test 
the former possibility and the detection of high-z massive 



clusters has bee n used to constrain the level of primordial 
non-gaussianity (Jimenez & Verde 2009; Holz & PerlmuttM| 
JOlO; Cayon et al. 2010; Grossi et al. 2007; Sartoris et al] 
[2OIOI : iHovle et al.l |2010|). In this paper, instead, we inves- 
tigate the second option: we point out that the presence of 
high-z massive clusters can be used to test the cosmological 
constant scenario as compared to dynamical coupled dark 
energy models. We therefore keep gaussian initial conditions 
and instead allow for the possibility that dark energy is not 
a cosmological constant but rather a dynamical scalar field 
that can interact with cold dark matter (CDM). 

We carry out N-Body simulations for a set of models of 
interacting dark energy with either constant or time depen- 
dent couplings. We then compute the cumulative halo mass 
functions of the structures that form within these theories 
at different redshifts, and compare them with the outcomes 
of the standard ACDM model. In this way we are able to 
show that the probability of finding massive clusters at high 
redshifts is significantly enhanced in coupled dark energy 
cosmologies with respect to ACDM. 



2 COUPLED DARK ENERGY 

We consider coupled dark energy models where t he acceler- 
ation of the Universe is driv en by a scalar field (j) l|Wetterichl 
ll988l : lRatra fc Peebleslll988h interacting with the CDM fluid 
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according to the equations: 



3 SIMULATIONS 



p'^ + ZHpcf, = +I3(j}'pc 



(1) 
(2) 



with pc the energy density of CDM and where a prime de- 
notes a derivative with respect to conformal time. The di- 
mensionless function 



/3: 



dlnrric 
d(l) 



(3) 



where rric is the mass of a CDM particle, fully specifies the 
interaction. Note that we define densities as p = SvrGp and 
that the scalar field (j> is normalized in units of the reduced 
Planck mass M = (87rG)-^''^ The scalar field energy density 
is p0 = §^ + U{(j)) and we consider two choices of potentials: 



C/(0) 



Uoe 



(4) 
(5) 



where Uq and a are constants. In the present work we limit 
our investigation to the models specified in Table [1] one 
model with a potential given by Q and with a constant 
coupling P — 0.2, labelled as "RP", and two models with 
an exponential potential ((5]) with constant coupling /3 = 0.2 
and variable coupHng /?(<;/>) = 0.4 ■ e^*^, labell ed "EXP005" 
and "EXP010a2", respectively fsee lBaldilbOlOl . for a detailed 
discussion on variable coupling models). 

Coupled dark energy models have b e en widely 
studi e d in the literature JWetterich 19951: Amendolal 
I2OOOI : iMangano et al. I I2OO3I : iPettorino fc Baccigalupil |2008| . 
and references therein). The interaction imprints spe- 
cific features on the growth of cosmic structures, as 
it h as been shown both wi t hin s p herical collapse mod- 
els (Wintergerst fc Pettorinol I2OI0I : iMainini fc Bonomettol 
|200(: ) and at the non linear level by means of N-body sim - 
ulations for constant (Baldi ct al. 2010: Maccio et al.ll2004l ') 
and variable couplings I3{(f>) ( Baldi, 2010, '). In particular, the 
interaction determines an enhancement of structure forma- 
tion due to the presence of a long range fifth-force acting 
between coupled massive particles. In addition to this effect, 
conservation of momentum turns into an extra acceleration 
directly proportional to the peculiar velocity of coupled par- 
ticles; at the nonlinear level, it has been shown that this 
velocity dependent term can determine a reduction of the 
concentration of halos if the coupling functio n I3{(fi) is posi- 
tive (see again iBaldi et al.ll2010l : lBaldill2010l '). These effects 
are encoded in the modified newtonian acceleration equation 
for coupled dark matter particles: 



"•"Z^ |^f,,|3 



(6) 



where i and j are indices that span over all the particles of 
the simulation, p is the momentum of the i-th particle, and 
Gij is the effective gravitational constant between the i-th 
and the j-th coupled particles 

G., =Giv [1 + 2/3^(0)], (7) 
where Gjv is the usual Newtonian value. 



In this Letter, we estimate how the probability of finding 
very massive clus ters comparable to t he XMMU J 2 235.3- 
2557 detection bv lJee et al.l (|2009l ) and lRosati et all (|2009D 
at 2 ~ 1.4 is modified in the presence of a coupling be- 
tween dark energy and CDM as compared to the standard 
ACDM case. With this aim, we use the mo dified version o f 
the cosmological N-body code GADGET-2 llSpringelllioOsh 
presented in lBaldi et all l|2010l ') and iBaldil (|2010l '). to which 
we refer for further details, to investigate nonlinear struc- 
ture formation within the models described in Table [1] We 
make use of the following sets of simulations: 

• two of t he high-reso l ution hydrodynamical simulations 
presented in lBaldi et al.l (|2010l ) for the case of a scalar field 
with an inverse-power law self interaction potential with 
a constant coupling /3 = 0.2 and for a reference ACDM 
model, both with cos mological parameters in accordance 
with WMAP5 results l|Komatsu et al.ll2009l '): 

• three new simulations at lower resolution in a larger 
cosmological box for the case of an exponential self inter- 
action potential with either constant (/? = 0.2) or variable 
(/3 — 0.4 • e^*^) coupling and for a reference ACDM cosmol- 
ogy, all with cos mological param eters in accordance with 
WMAP7 results (iKomatsu et al.| [2010'l. 

The former set of simulations includes the cosmological frac- 
tion of uncoupled baryons on which hydrodynamical forces 
are computed with the Smoothed Particl e Hydrodynamics 
l|Springel fc Hernauistl2002l : ISpringeill2005l ) algorithm, while 
the latter ones are pure CDM simulations. 

The two sets of simulations start at 2 = 60 with ex- 
actly the same initial conditions for each set of runs. This 
is a conservative setup, since the effects of enhanced growth 
that we are investigating would be more pronounced if a 
normalization at decoupling (2 ~ 1100) had been adopted. 
We remark that this normalization is different f rom the one 
ado£ted in most of the simulations discussed in lBaldi et al.l 
(j2O10l ) and lBaldil l|2010l ): in these works the amplitude of Un- 
ear density fluctuations was normalized in order to have the 
same erg at the present time. In the present study, instead, 
we choose all the models with the same initial normaliza- 
tion of the power spectrum, which will necessarily result in 
different values of erg at 2 = for the different cosmological 
models. It is important to notice that all our models start at 
high redshift with a normalization of the Power Spectrum 
which is in full accordance with the latest WMAP7 deter- 
mination of the scalar perturbat ions amplitude from CMB 
data alone (iKomatsu et al.ll2010l ). 

Given this setup, we study the halo mass function at 
different redshifts for the groups identified in each simula- 
tion with a Friends-of-Friends (FoF) algorithm with a linking 
length A — 0.2 x d, where d is the mean particle spacing. We 
show that the interaction between dark energy and CDM 
results in a significantly larger number of massive halos at 
any epoch, with respect to ACDM. 



4 RESULTS 

In Fig. ([1} we show the evolution of the halo mass function 
at different redshifts in the three low resolution simulations 
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Model 


U((l>) 


ex 


R 

M 


Box Size 
(h-i Mpc) 


Number 
of particles 


Mi 

{H-^Mq) 




(/i~i kpc) 


ACDM (low) 









320 


256=* 




1.47 X 10" 


25.0 


ACDM (high) 









80 


2 X 512^ 


4.7 X lO'^ 


2.3 X 10** 


3.5 


RP 


— 


0.143 


0.2 


320 


2 X 512^ 


4.7 X lO'^ 


2.3 X 10** 


3.5 


EXP005 




0.1 


0.2 


80 


256^ 




1.47 X 10" 


25.0 


EXP010a2 




0.1 


0.4 ■ e^'t' 


320 


256^ 




1.47 X 10" 


25.0 



Table 1. List of the different simulations discussed in the present work, performed with our modified version of GADGET-2. The low 
resolution simulations include only CDM particles while the two high resolution runs include uncoupled baryons with hydrodynamical 
forces. The last column of the table reports the gravitational softening ts used in each simulation. 




Figure 1. Cumulative mass functions for the low-resolution simulations of the coupled DE models with constant (red) and variable 
(green) couplings, and for the standard ACDM cosmology (black). The four plots correspond to different redshifts and show in the 
bottom panel the ratio of the halo number density to the ACDM case. The mass functions are based on the groups identified with a FoF 
algorithm and the masses quoted in the plots are FoF masses. Each plot also reports the mass of the most massive halo found at a given 
redshift in each simulation, which clearly shows how more massive structures are expected to form at any cosmological epoch in coupled 
dark energy models as compared with ACDM. 



described in Table [T] The bottom panel of each plot shows 
the enhancement of halo number density in coupled dark 
energy models with respect to ACDM. As it clearly appears 
in all the plots, the number density of halos of any mass 
is larger in coupled dark energy models, both for constant 
and variable couplings, as compared to ACDM, and the ef- 



fect sign ificantly grows w i th ma ss (the decreasing behavior 
found in Manera fc Motal ||2006| ) can be explained in view of 
IWintergerst fc Pettorinol l|2010l )). In particular, at 2 ~ 1.5 
the halo number density at the high-mass end of our simu- 
lated mass functions exceeds the ACDM value by a factor 
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Figure 2. Cumulative mass function at z = 1.5 for the coupled 
DE model RP of Table [l] (orange) compared to ACDM (black) 
from a high-resolution hydrodynamical simulation in a box of 
80 comoving Mpc. The enhancement of the mass function in 
this set of simulations is slightly lower than in pure CDM sim- 
ulations due to the presence of a fraction of uncoupled baryons. 
The bottom panel shows the enhancement of halo number density 
relative to ACDM. 

~ 10 and ~ 3 for constant and variable couplings, respec- 
tively. 

In Fig. ([2]) we plot the same quantities for our hydro- 
dynamical high-resolution simulations, which confirm the 
trend found in the larger simulation box, although the en- 
hancement is shghtly weaker in this case due to the inclusion 
of the fraction of uncoupled baryons. 

For all simulations, the most massive halo forming in 
the simulated box has a higher mass for coupled dark energy 
cosmologies than for ACDM. In this respect, it is important 
to stress once more here that all the simulations start from 
the very same initial conditions, and therefore this effect is 
due only to the different physics induced by the coupling and 
is not affected by statistical differences in the random real- 
ization of the initial density fields. Due to our normalization 
choice, the different models will also be characterized by dif- 
ferent values of ag at z = 0. In order to quantify this effect 
we have computed 0-8(2: = 0) from our two hydrodynamical 
simulations. The maximum increase in as is obtained, as ex- 
pected, for the highest value of the coupling: (Tg(/3 — 0.2) ~ 
0.91. Compared to the ACDM value crg(ACDM) ~ 0.76, 
this corresponds to a maximum relative increase of roughly 
20%, still in reasonable agreement wit h present available 
measurements of gg at z = (see e.g. iMantz et al.l l2009l : 
IVikhlinin eral]|2009l ). We stress again that different values 
of (Tg at the present time arise from the same normaliza- 
tion at high redshift. In other words, different couplings are 
able to evolve the same initial normalization into different 
(Tg today. As a consequence, a mismatch between the value 
of the perturbations amplitude at CMB and local measure- 
ments of as based solely on low redshift probes, could pro- 
vide a further possible observational test for coupled dark 
energy models, breaking degeneracies with ACDM cosmolo- 
gies. Interestingly, some mismatch between the value of erg 
extrapolated from CMB data under the assumption of a 
standard ACDM cosmology and p ure local measu r ement s 
has been recently claimed by e.g. lFeldman et ahl l|2003h : 



iReichardt eFal] (|2009l '): IWatkins et all l|2009l ). although still 
with a low statistical significance. 

We note also that the effect is expected to be model- 
dependent and possibly weaker for other viable models as 
e.g. variable coupling scenarios. 



5 CONCLUSIONS 

We conclude that coupled dark energy can significantly in- 
crease the probability to detect massive clusters at high red- 
shift with respect to a ACDM model. The effect shows a 
quite strong dependence on the halo mass and is found to 
be larger for higher halo masses. Note that similar coupled 
dark energy models may also predict the existence of struc- 
tures at very large scales jWintergerst et al]|2010l ). 

For the case of a constant coupling we have assumed 
the largest possible value of the coupling which is allowed 
according to present constraints for this class of mode ls 
iBean et al.ll2008l : iLa Vacca et"alll2009l : iBaldi fc Vieill2010l ') . 
In this case the number density of halos at z ~ 1.5 is found 
to be at least a factor ~ 10 higher than in ACDM cosmolo- 
gies for halos of masses larger than M « 10^* Mq. This 
is a conservative estimate and the enhancement could in- 
crease in case the models were normalized further in the 
past (e.g. at decoupling) rather than at 2: ~ 60 as considered 
in the present simulations. 

We have also presented results for one choice of variable 
coupling models. Despite the large value of the interaction 
strength at recent cosmological epochs, such models have a 
weaker impact on the background expansion history as com- 
pared to constant coupling models and are therefore easier 
to reconcile with present cosmological constraints. For the 
specific choice of variable coupling that we have considered 
in this work, which has recently been shown to also mitigate 
other possible tensions between the ACDM cosmology and 
observations at small scales (as e.g. the cluster baryon frac - 
tion or the so called "cusp-core" problem, see lBaldll2010l ). 
the total effect on the number counts is found to be some- 
what weaker than in the constant coupling scenario, due to 
the low value of the coupling at early times. Nevertheless, 
also in this model the probability of finding massive clusters 
at z ^ 1.4 is enhanced by at least a factor 2 — 3 with respect 
to ACDM. 

These estimates are based on the ratio of the simulated 
halo mass functions at their high-mass end in the different 
models. Although the most massive groups found in our 
simulations at z ~ 1.4 do not reach the m ass of the high-z 
clust er XMMU J2 2 35. 3-2 557 identified bv lJee et al.1 (|2009l ) 
and iRosati et all (|2009l '). the extrapolation to slightly 
higher masses of the enhancement found in the number 
density at the high-mass end of our halo catalogs is a 
conservative approach since the effect is found to be quite 
strongly increasing with mass (see Figs. Ill2p . It is also 
interesting to notice that the mass of the most massive 
halo forming in coupled dark energy cosmologies, both for 
constant and variable couplings, is systematically larger 
than for the corresponding ACDM cosmology. The related 
increase of as does not exceed 20% for the largest possible 
coupling and could be significantly reduced for other viable 
scenarios or by performing an optimization of the full set 
of cosmological parameters. In any case, such an increase 
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concerns only the value of as at low redshifts, while all 
the models are consistent at high redshift with the present 
CMB data. One may want to compare the effects found 
in the present analysis for coupled dark energy with those 
of a ACDM model with a higher ug: wc have noted that 
such a degeneracy may be broken by directly measuring the 
redshift evolution of as, thereby providing an additional 
observational test for coupled dark energy as compared to 
ACDM. 

We have therefore shown that both constant and vari- 
able coupling models, also for different scalar potentials, en- 
hance at any cosmological epoch the cumulative halo mass 
function and consequently the probability to detect halos 
of any given mass in volume limited surveys with respect to 
the standard ACDM cosmology. Future detection of massive 
clusters at high redshift can therefore be used to disentan- 
gle a cosmological constant from dynamical coupled dark 
energy models. 
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